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Azimuth multipoles vs minimum-bias jets in 2D angular correlations on r\ and 
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Angular correlations measured in p-p and heavy ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) include a same-side (SS) 2D peak. In peripheral A-A and p-p collisions the SS peak 
properties are consistent with predicted minimum-bias jet correlations. However, in more-central 
Au-Au collisions the SS peak becomes elongated on pseudorapidity n. Arguments have been pro- 
posed to explain the SS peak n elongation in terms of possibly-fluctuating initial-state geometry 
multipoles coupled with radial flow to produce final-state momentum-space multipoles. Such argu- 
ments are based on Fourier decomposition of 2D angular correlations projected onto ID azimuth. In 
this analysis we show that measured correlation structure on rj (large curvatures) establishes a clear 
distinction between the SS 2D (jet) peak and ID multipoles. Measured 2D peak systematics can 
predict inferred ID Fourier amplitudes interpreted as "higher harmonic flows." But ID Fourier am- 
plitudes alone cannot describe 2D angular correlations. The SS 2D peak remains a unique structure 
which can be interpreted in terms of parton scattering and fragmentation in all cases. 

PACS numbers: 12.38.Qk, 13.87.Fh, 25.75.Ag, 25.75.Bh, 25.75.Ld, 25.75.Nq 



I. INTRODUCTION 

Data from nuclear collisions at the Relativistic Heavy 
Ion Collider (RHIC) have been interpreted within a hy- 
drodynamic (hydro) context to demonstrate the forma- 
tion of a thermalized, flowing partonic medium with 
small viscosity 0, 0|- However, alternative analysis of 
spectrum and correlation data reveals a minimum-bias 
jet or minijet contribution whose variation with A-A cen- 
trality and collision energy seems to conflict with hydro 
expectations 043 ■ A large-amplitude 2D peak at the ori- 
gin in angular correlations on pseudorapidity r\ and az- 
imuth (f> expected as a jet manifestation persists even in 
central Au-Au collisions, albeit the 2D peak is elongated 
on rj relative to a nominally symmetric jet cone 

The competition between minimum-bias jets and flows 
to explain the final-state structure of RHIC collisions 
has recently intensified and now focuses on the SS 2D 
peak in p t -integral angular correlations. The minijet in- 
terpretation of the SS peak in more-central A-A colli- 
sions 0, 3] has been challenged. The peak has been 
reinterpreted in terms of "triangular flow" resulting from 
conjectured initial-state transverse geometry fluctuations 
(i.e., sextupole and other azimuth multipoles) coupled 
to radial expansion or as resulting from interac- 

tion of initial-state "glasma flux tubes" coupled to radial 
flow [ll|, [HI • Interpretations based on initial-state A-A 
overlap geometry (the subject of this article) typically in- 
voke a Fourier series to describe 2D angular correlations 
projected onto ID azimuth. The Fourier sinusoids are in- 
terpreted to represent "higher harmonic flows" [Io|, HH ■ 

However, the ID projection onto azimuth which is cen- 
tral to such arguments discards critical information con- 
tained in the pseudorapidity structure. A discrete Fourier 
series can describe any structure projected onto periodic 
azimuth. A Gaussian peak narrow on azimuth must con- 
tribute significantly to several ID Fourier terms. But 
the SS 2D peak contribution should (and can) be dis- 



tinguished from other 2D structure based on its strong 
pseudorapidity variation (curvature). 

In this analysis we invoke measured systematics of 
the minimum-bias (pt-integral) SS 2D peak for 200 GeV 
Au-Au collisions to demonstrate that recently-reported 
v m {2} measurements correspond to properties of the SS 
peak. From systematics of the peak properties we can 
predict any higher harmonic flow v m {2} (Fourier coef- 
ficient) as a function of centrality and rj exclusion cuts 
meant to remove "nonflow" from v m measurements. 

In p-p and more-peripheral A-A collisions the SS 2D 
peak conforms closely to expectations for minimum-bias 
jets [1, [T3 - [T7j . Based on comparison of spectra, cor- 
relations and pQCD calculations we conclude that the 
SS peak continues to represent jet production even in 
more-central Au-Au collisions 0t3]- The most likely in- 
terpretation of the SS peak mechanism in more-central 
A-A collisions remains parton scattering and fragmenta- 
tion, with modification of fragmentation including r/ elon- 
gation (polarization) of lower-momentum jet fragments. 
Thus, jet production is the mechanism behind the mul- 
tipoles recently interpreted as higher harmonic flows. 

This article is arranged as follows: We review analysis 
methods and possible correlation mechanisms in Sees. [TT1 
through llVI We summarize measured 2D angular correla- 
tions from 200 GeV Au-Au collisions and the systematics 
of minijet structure and the nonjet azimuth quadrupole 
in Sec. |V] We discuss Fourier series analysis and possible 
confusion between minijets and a nonjet quadrupole aris- 
ing from some analysis methods in Sees. IVII and lVIIl We 
then review conjectured initial-state geometry structure 
and related measures (triangular flow, higher harmonics) 
in Sees. I Villi and IIXI In Sec. [X] we present quantitative 
relations between minijet structure and azimuth multi- 
poles inferred from various analysis techniques, with a 
direct comparison between higher multipoles predicted 
from 200 GeV minijet systematics and a recent measure- 
ment of higher harmonic flows from the LHC. 
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II. ANALYSIS METHODS 

We briefly introduce correlation analysis methods ap- 
plied to nuclear collisions at the RHIC. Method details 
are described in Refs. |E II El EE EESI1 • Topics include 
A-A collision geometry, correlation measures and 2D cor- 
relation spaces. A detailed discussion of initial-state A-A 
geometry is presented in Sec. IIV1 a two-component an- 
gular correlations data model is presented in Sec. [V] and 
Fourier series analysis relevant to ID azimuth correla- 
tions is reviewed in Sec. IVI1 



A. Initial-state A-A geometry 

Initial-state (IS) A-A collision geometry is described 
by the Glauber model relating the A-A differential cross 
section to participant-nucleon number N part and N-N 
binary-collision number A^ in 22] . A derived participant- 
nucleon mean path length v = 2N{,i n /N part can also be 
defined. Through the measured A-A differential cross 
section on charged-hadron multiplicity n c h within some 
angular acceptance the Glauber model parameters are 
related to observed n c h- 

Optical e opt and Monte Carlo e^/c (HI eccentric- 
ities have been invoked to model the IS A-A overlap 
eccentricity required for interpretation of the final-state 
(FS) azimuth quadrupole measured by x>2- (-opt assumes 
a smooth matter distribution across nuclei whereas cmc 
assumes that point-like participant nucleons are the de- 
termining elements. A priori support for e opt derives 
from a conjecture that the nonjet azimuth quadrupole 
emerges from interactions at small x < 0.01 where one 
might expect onset of a smooth, saturated glue system 
(e.g. Glasma) |24|. A posteriori support for e op t is sug- 
gested by a simple systematic trend observed for the non- 
jet quadrupole fSec lVCj) . 



B. Two-particle correlation measures 

Two-particle correlations are structures in pair- 
density distributions on six-dimensional momentum 
space (pti ,rji,(pi, pt2 , f]i , <f>2 ) ■ We visualize correlation 
structure in 2D subspaces (pt,Pt) and (t]a,<Pa) (defined 
below) which retain almost all structure within a lim- 
ited r\ acceptance such as the STAR Time Projection 
Chamber (TPC). Correlations can be measured with per- 
particle statistic Ap/^/p re f — po ((r) — 1), where Ap = 
Psib — Pref is the correlated-pair density, p s n, is the sibling 
(same-event) pair density, p re f is the reference- or mixed- 
pair density, (r) is the mean sibling/mixed pair-number 
ratio, and prefactor po = n c h/ Ar] A(f> sa d 2 n c h/ drjd(f> is 
the charged-particle 2D angular density averaged over 
angular acceptance ( Ary, A<j>) [E [25[ . Pair ratio r is av- 
eraged over kinematic bins (e.g. multiplicity, p tl vertex 
position). Factorization p re f sa p\ is assumed. 



The per-particle quadrupole component of 2D angu- 
lar correlations resulting from some analysis method is 
Aq {method} defined in Eq. ([5]). Correlations can also 
be measured with per-pair statistic Ap/p re f, including 
total azimuth quadrupole component f|{2} and higher 
multipoles v^ n {2}. Variation of per-pair correlation mea- 
sures with A-A centrality is typically dominated by a 
trivial l/n c h trend, or in the case of v m a l/y/n c h trend. 

C. Number correlations on (pt,Pt) or (yt,yt) 

2D correlations on pt or transverse rapidity yt = 
hi[(pt + Jnt)/m 7r ] (m-jr is assumed for unidentified 
hadrons) are complementary to 4D angular correlations 
in 6D two-particle momentum space. Correlations on an- 
gle differences (t?a,^a) and momenta (pt,Pt) [or rapid- 
ity (yt,yt)] can be defined for like-sign (LS) and unlikc- 
sign (US) charge combinations and also for same-side 
(SS) and away-side (AS) azimuth subregions of angu- 
lar correlations (defined below). Manifestations of differ- 
ent correlation mechanisms (e.g. so-called soft and hard 
components, Sec. IIIip can be clearly distinguished in the 
four combinations of charge-pair type and azimuth sub- 
space, with distinctive forms for each of the LS and US 
charge combinations and for SS and AS azimuth sub- 
spaces [bH EB] • Any conjectured correlation mechanism 
must accommodate all such observed systematic trends. 

D. Number and pt angular correlations on (rjA,<t>A.) 

Angular correlations can be formed by integrating over 
the entire (pt,Pt) pair acceptance (minimum-bias angular 
correlations) or over subregions [lj, EH • Examples of the 
latter include "trigger-associated" dihadron correlations 
resulting from asymmetric cuts on (pt,Pt) [261 ] - 

Two-particle angular correlations are defined on 4D 
momentum subspace (771, 772, <f>i, <f>2)- Within acceptance 
intervals where correlation structure is approximately in- 
variant on mean polar or azimuth angle (e.g. rys = 771 +772) 
angular correlations can be projected by averaging onto 
difference variables (e.g. t/a = rji — rj 2 ) without loss of in- 
formation to form angular autocorrelations [IE [25| ■ 2D 
subspace (r?A, <M) 1S then visualized. Symbol Aa; denotes 
the detector acceptance on parameter x. The pair angu- 
lar acceptance on azimuth can be separated into a same- 
side (SS) region (|0a| < tt/2) and an away-side (AS) 
region (\4>a\ > tt/2). The SS region includes mirojet cor- 
relations (hadron pairs within single jets) , while the AS 
region includes internet correlations (hadron pairs from 
back-to-back jet pairs). 

Attempts to isolate 2D or ID correlation structure 
from a (large) combinatoric background generally fol- 
low one of two methods: (a) model fits to 2D his- 
tograms [E H and (b) ZYAM subtraction from ID di- 
hadron correlations on azimuth [27j |. 2D angular corre- 
lations on (t/Aj^a) are observed to include a few ele- 
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merits accurately described by simple functional forms, 
as described in Sec [V] Two of the elements have been 
interpreted in terms of jet correlations 0, H| • 



III. HADRON PRODUCTION MODELS 

Several classes of hadron production models are in- 
voked to describe high-energy nuclear collisions at the 
RHIC and LHC. One class consists of N-N superposition 
models, the limiting case being Glauber linear superposi- 
tion (GLS). A second class is based on the limiting case of 
a homogeneous bulk medium produced through substan- 
tial parton and/or hadron rescattering with significant 
equilibration of low-energy partons. High-energy par- 
tons are expected to probe thermalized-medium proper- 
ties via jet structure modification. A third class incorpo- 
rates "nonsmooth" initial conditions plus hydrodynamic 
evolution to produce novel correlation structure in the 
final state. 



A. N-N superposition models 



angular correlations is an azimuth quadrupole described 
by function cos(2</>a) and interpreted as elliptic flow [35| . 

Medium formation and thermalization rely on 
hadron [U IH[ and possible parton rescattering. 
Hadrons emerge from the bulk medium through a freeze- 
out process exhibiting radial [3l|,|36| as well as elliptic (35j 
flow. Hadrochemical trends and pt spectrum structure 
are in terp reted to indicate a thermal system supporting 

flows jUt HU. 



C. Initial-state geometry plus radial flow 

More recently, higher multipoles in initial-state A-A 
geometry coupled to radial flow in a bulk medium have 
been proposed as a possible mechanism for some final- 
state correlation structure. For instance, the SS 2D peak 
is conjectured to be a consequence of initial-state geom- 
etry multipoles [9(, geometry fluctuations [8j or Glasma 



flux tubes [12|. Theoretical models are represente d by 



Monte Carlos such as AMPT Q and NexSpheRIO [37|. 
IS geometry models compete directly with N-N superpo- 
sition models. That dichotomy is a subject of this article. 



Hadron production from N-N collisions, as in unmodi- 
fied single collisions (linear) , or collisions modified by the 
A-A environment (nonlinear), is superposed according to 
the Glauber model to describe hadron production in the 
A-A final state. Hadron production in N-N collisions is 
described by a theoretical two-component model includ- 
ing longitudinal projectile fragmentation and transverse 
scattered parton fragmentation [28U30j . Modified hadron 
production in A-A collisions [5| can be described by al- 
teration of the pQCD fragmentation (parton splitting) 
process 0]. 

Longitudinal nucleoli fragmentation results from "soft" 
four-momentum transfers between projectiles. Trans- 
verse parton fragmentation results from "hard" momen- 
tum transfers between constituent partons at some mo- 
mentum fractions x\, X2 leading to dijet production. The 
terms "soft" and "hard" refer to the IS four-momentum 
transfer, not to the momenta of FS hadrons. Correlations 
from the soft component play a negligible role in more- 
central A-A collisions, but the soft component remains 
the dominant single-particle hadron production mecha- 
nism. The nonjet azimuth quadrupole is introduced in 
Sec. [V] as a "third component," first observed in more- 
central A-A collisions. 



B. Bulk- medium formation and freezeout 

A second model class emerging from lower-energy 
heavy ion programs [3lj describes production of a ther- 
malized bulk medium (hadronic and/or partonic) in 
which the dominant collision mechanism is hydrody- 
namic flows 3J|. The signature manifestation in A-A 



IV. INITIAL-STATE A-A GEOMETRY 

Initial-state A-A overlap geometry was first modeled 
by smooth nuclear-matter densities to obtain "optical" 
eccentricity e opt as a function of centrality. Motivated 
in part by interest in alternative interpretations of an- 
gular correlation features possibly attributed to jets the 
Glauber Monte Carlo was extended to derive tMC from 
the distribution of point- like participant nucleons [23| . 
It is conjectured that IS geometry fluctuations coupled 
with radial flow may produce flow structure appearing as 
jet-like FS correlations The role of conjectured fluc- 
tuations in IS phase-space geometry may be relevant to 
interpretation of such jet-like structure. 



A. IS transverse phase space 

IS structure is present in both momentum space and 
configuration space, whereas FS correlation structure is 
observed only in momentum space. IS momentum-space 
structure (e.g., scattered-parton distributions) may be 
transported with some fidelity to FS momentum space 
via a separately- measured process (e.g., parton fragmen- 
tation to jets) , whereas IS configuration-space structure 
(e.g., conjectured IS geometry fluctuations) requires cou- 
pling to flows for manifestation in FS momentum space. 

Minimum-bias parton scattering is a common feature 
of IS momentum space in all high-energy nuclear col- 
lisions. The scattered-parton distribution represents a 
large range of projectile-nucleon momentum fraction x, 
whereas the IS nonjet quadrupole may emerge only from 
low x < 0.01. Whether large-angle-scattered partons 
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near mid-rapidity fragment to detectable FS jets in A-A 
collisions or whether they thermalize rapidly to drive ra- 
dial flow is a major issue for RHIC collisions. 



eccentricity for 200 GeV Au-Au is parametrized by [2(| 



C2 ' opt ~ 5A 



log 
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3N, 
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log 
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/1136 



0.78 



(2) 



B. IS azimuth power spectrum 

The IS A-A transverse geometry (at x « 1) may have 
three components: a static ("optical") part at fixed im- 
pact parameter 6, a contribution from fluctuating b, both 
represented by even azimuth multipoles, and a conjec- 
tured stochastic part (point-like participant sampling) 
represented by a white-noise spectrum including even and 
odd multipoles. 

A-A initial-state azimuth structure modeled at I « 1 
by a Glauber Monte Carlo is described by a participant- 
nucleon autocorrelation [ijj]- For non-central A-A colli- 
sions the autocorrelation on azimuth difference 4>& in- 
cludes a few even-m sinusoids dominated by m = 2 
and phase-correlated with impact parameter b, a uni- 
form combinatoric background and a delta function oc 
Npart (self pairs) representing participant-nucleon sam- 
pling noise. By the Wiener-Khinchine theorem the 
Fourier transform of the azimuth autocorrelation is a 
power spectrum represented by eccentricity elements 



= N part e 2 ml with per- pair eccentricity measures \24 



c rn,MC 



SeL m odd. 



m even 



(1) 




FIG. 1: Centrality trends for optical and Monte Carlo 
Glauber initial-state azimuth multipoles with m = 2, 3, on 
participant-projectile-nucleon number (left panel) and binary 
N-N collisions v per participant-nucleon pair (right panel). 

By comparing the trends in Fig. Q] with published 
Monte Carlo results (e.g., Ref. Q - Fig. 2) we find 
that the stochastic contribution to the power spectrum 
for m = 2, 3 is consistent with O(l) x 2/N part , with 
0(1) ~ 0.5-2 as expected for a per-pair fluctuation mea- 
sure and Poisson statistics. Whether point-like sampling 
represents IS geometry with significant manifestations in 
FS correlation structure is an open question. 



Eccentricity e^opt represents the "elliptical" A-A over- 
lap region for fixed b, and cr\ represents the eccentricity 
variance due to event-wise fluctuations in b. Monte Carlo 
random sampling generates a power spectrum Se^ oc 
1/ Np ar t approximately uniform on m corresponding to 
the self-pair contribution « N part 5(4>A) hi the azimuth 
autocorrelation. For a stochastic process there should 
be no phase relation between noise amplitudes Se^ and 
impact parameter b. All higher m are present in the IS 
Monte Carlo spectrum and might appear in the final state 
to some extent if Monte Carlo sampling at i ~ 1 were 
a legitimate model of IS geometry relevant to FS hadron 
production for x < 0.01 and r\ = 0. 



C. Azimuth power spectrum centrality trends 



V. TWO-COMPONENT 2D DATA MODEL 

Correlation data in the form of 2D angular autocor- 
relations [l|| can be obtained for various pt cut con- 
figurations, includingso-called "trigger-associated" di- 
hadron correlations [26| and minimum-bias (p t -integral) 
correlations. For simplicity of illustration we restrict 
to p t -integral data from 200 GeV Au-Au collisions 
This exercise illustrates construction of a necessary and 
sufficient (N-S) 2D mathematical data model based en- 
tirely on data phenomenology (see App. [K§. The three 
model elements are subsequently interpreted physically 
in terms of longitudinal projectile fragmentation (soft 
component), transverse parton fragmentation (minijets) 
and a nonjet azimuth quadrupole whose physical origin 
is in question. 



Figure [T] shows centrality trends for m — 2, 3 IS ge- 
ometry power spectrum elements on participant-nucleon 
number N part (left panel) and mean participant path- 
length v (right panel). Plotted are optical eccentric- 
ity £2, opt (solid curve), Monte-Carlo eccentricity e2,MC 
(dash-dotted curve) and so-called "triangularity" 5e 3 
(dashed curve). From Eq. ([T]) (and ignoring a possi- 
ble of 2 contribution) we have e\ MC = e\ opt -I- 5e\ with 
5e% w 4/Npart and e\ MC = 5e§ w 4/N part . The optical 



A. Angular correlation model function 

2D data histograms on (t?a>^a) are n °t generally fac- 
torizable. However, minimum-bias histogram data ob- 
tained from RHIC p-p and A-A collisions can be repre- 
sented as the sum of a few factored terms 

h(rjA,(pA) = ^2 /«( 7 ? A ) 9n(M- (3) 

n 
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For minimum-bias data the series accurately represents 
all information in the data histogram with a few sim- 
ple model functions. For some terms one factor may be 
approximately constant, further simplifying the model. 

Minimum-bias data from Au-Au collisions are de- 
scribed by three main elements: (a) a same-side (SS) 2D 
peak at the origin on (j7Aj<^a) well approximated by a 
2D Gaussian for all minimum-bias data, (b) an away-side 
(AS) ID peak on azimuth or "ridge" well approximated 
by AS azimuth dipolc [1 — cos(</>a)]/2 for all minimum- 
bias data and uniform to a few percent on ?7a (hav- 
ing negligible curvature) , and (c) an azimuth quadrupole 
cos(2(/>a) also uniform on r]A to a few percent over the 
full angular acceptance of the STAR TPC. Model el- 
ements (a) and (b) together have been interpreted as 
minimum-bias jets or "minijets" 0- Element (c), de- 
scribed as the (nonjet) azimuth quadrupole, is conven- 
tionally attributed to elliptic flow, a hydrodynamic phe- 
nomenon |35j (but see App. IB"]) . 

The 2D model function [Eq. equivalent] applicable 
to more-central A-A collisions is 0, 0, H3] 



^Pref 



= A + A 2 d exp 



VA 



+ ^1d{1 + cos(0a -tt)}/2+/1q2cos(2 0a).(4) 



A ID Gaussian on t^a modeling projectile nucleon frag- 
mentation (soft component, negligible in more-central 
Au-Au collisions) and a 2D exponential modeling quan- 
tum correlations (HBT) and electron pairs (extremely 
narrow in more-central Au-Au collisions) are omitted for 
simplicity in discussion of more-central A-A collisions. 
Quadrupole measure Aq is statistically consistent with 
jet measures A 2 d and Ajj, permitting quantitative com- 
parisons between jet and nonjet-quadrupole systematics. 

Figure [5] shows examples of 2D angular correlations 
from four centralities of 200 GeV Au-Au collisions based 
on fit parameters of Ref . |4J . The centrality values v cor- 
respond to (a) 1.25 (« N-N collisions), (b) 2.5, (c) 4.5 
and (d) 6 (b = 0). The histograms are plotted within 
the STAR TPC angular acceptance |?7a| < 2 convention- 
ally adopted for 2D correlation analysis. 2D correlation 
structure is accurately described by a simple mathemat- 
ical model, and 2D fit residuals (r.m.s. amplitude) are 
typically less than 1% of the SS 2D peak amplitude for 
more-central Au-Au collisions. Residuals for ID projec- 
tions onto azimuth are substantially smaller because of 
averaging. We now consider the details of minijet and 
nonjet quadrupole systematics in turn. 



B. Final-state minijet systematics 

Figure [3] summarizes preliminary fitted SS 2D and 
AS ID peak parameters vs centrality measure v within 
the nominal STAR TPC angular acceptance (Ary, A(f>) = 
(2, 2ir) A 2 d in the left panel is the fitted amplitude of 
the SS 2D Gaussian function. Its two r.m.s. peak widths 




FIG. 2: (Color online) Angular correlation histograms for four 
centralities of 200 GeV Au-Au collisions based on fit param- 
eters from Ref. [J]. The centralities measured by v are given 
by (a) 1.25 (~ N-N), (b) 2.5, (c) 4.5 and (d) 6 (6 = 0). 



are shown in the right panel. AS dipole amplitude Ad 
closely follows the SS 2D peak amplitude as expected 
for back-to-back dijets. There is smooth variation with 
centrality, but a "sharp transition" in SS 2D peak prop- 
erties occurs near centrality v w 3. Although the SS 2D 
peak becomes broad on 77A the peak curvature on 7?a re- 
mains large in all cases. Large curvature and unique pt 
correlation structure differentiate the SS 2D peak from 
conjectured flow mechanisms. 




FIG. 3: Left: Amplitude of the same-side 2D Gaussian A2D 
and away-side dipole Ad fitted to minimum-bias 2D angular 
correlation data from 200 GeV Au-Au collisions Right: 
Fitted peak widths for the same-side 2D Gaussian. GLS in- 
dicates a Glauber linear superposition reference extrapolated 
from measured p-p collisions [lq ]. 

The correlation parameters in Fig. [3] were combined 
with a pQCD parton spectrum to predict equivalent jet 
fragment yields as hadron spectrum hard components [6j . 
The predictions agree closely with previously-extracted 
spectrum hard components [5| . In turn, differential spec- 
trum hard components are described quantitatively by 
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a full pQCD calculation The agreement among 

pQCD theory, single-particle spectra and two-particle 
correlations within a common minijet framework pro- 
vides strong support for a minijet interpretation of the 
SS 2D peak for all Au-Au centralities. 



C. Final-state nonjet quadrupole systematics 

Figure [4] (left panel) shows centrality trends (solid 
curves) inferred from quadrupole data for 62 and 200 
GeV Au-Au collisions reported in Ref. [2(|. Per-particle 
quadrupole amplitude Aq is related to conventional mea- 
sure V2 by 



Agjmethod} = poifyvUmethod} 



(5) 



for V2 methods {2} w {EP} (nongraphical numerical 
methods) and {2D} (model fits to 2D histograms), where 
single-particle density po(b) is defined in Sec. Ill Bl Some 
of the 200 GeV Aq{2} data (open squares) fall above 
the upper plot boundary. v 2 {EP} data from 17 GeV 
Pb-Pb collisions are also included (solid triangles) [39| . 
The same centrality trend (solid curves) describes the 
Aq{2D} data over a large range of energies. 
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FIG. 4: Left: Quadrupole amplitudes for different collision 
systems and analysis methods. Solid curves are Aq{2D} from 
Ref. [2Cjl . Open squares are transformed «2{2} data from 
Ref. [231 . Solid triangles are similarly derived from V2{EP} 
data reported in Ref. [3IJ. The curves are described in the 
text. Right: Quadrupole data reported in the form of conven- 
tional elliptic flow measure v%. The curves are appropriately 
transformed from the straight lines in Fig. [5] (left panel). 

Figure @] (right panel) shows the energy systematics of 
quadrupole data in the form of conventional measure V2 , 
the square root of a per-pair correlation measure which 
tends to overemphasize small data values from systems 
with small particle multiplicities (e.g., lower collision en- 
ergies). The open circles are derived from published 
V2{EP} data, the solid points from Ref. [20j. Centrality 
choice b/bo w 0.5 (y ~ 4.3) minimizes the relative jet 
contribution to V2{EP}. The solid and dashed curves 
are derived from the straight lines in Fig. [5] (left panel) 
suitably transformed. 

Figure [5] (left panel) shows the energy systematics of 
quadrupole data in the form of per-particle measure Aq . 
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FIG. 5: Left: Collision energy dependence of quadrupole am- 
plitudes Aq for V2 data from Ref. [HI (open circles) and 
Aq{2D} reported in Ref. 20]. The solid line is from Eq. ©. 
The dashed line describes V2 data in Aq format at lower ener- 
gies [? ]. Right: Universal systematics for nonjet quadrupole 
Aq{2D} reported in Ref. [201 ] (bold solid line) compared to 
Aq{2} and Aq{EP} data presented in Fig. [4] (left panel). 



The Aq trend above 13 GeV is described by the solid 
line ^q(v^vtv; b/b « 4.5) = 0.13 R(^sWJt) 0, with 



R(V^) - In (ViWl3.5 GeV) . 



(6) 



The dashed line describing data below 13 GeV is 
0.008 In (y/J/vTv/3.2 GeV) describing the well-known 
transition (sign change) from squeezeout (due to partici- 
pant shadowing) to in-plane expansion. That panel sug- 
gests qualitatively different quadrupole production mech- 
anisms below and above 13 GeV, in contrast to Fig. 0] 
(right panel) which might suggest continuation of Be- 
valac/AGS projectile- nucleon (nucleon cluster) collectiv- 
ity to RHIC and LHC energies. 

Figure [5] (right panel) summarizes the measured cen- 
trality and energy dependence of nonjet quadrupole am- 
plitude Aq{2D}. All p t -integral nonjet quadrupole data 
from Au-Au collisions are accurately summarized above 
13 GeV by [H 



A Q {2D} = 0.0045i?(V5^) N bin e 



2, opt 



(7) 



defining the bold solid line which transforms to the solid 
and dashed curves in Fig. [4] (left panel). The simple lin- 



not e MC (dash-dotted curve). 



ear trend applies to e 2 opt 

The centrality trend of ^Q{2D}'on b/b (Fig. [4] - left 
panel) is Gaussian to good approximation, independent 
of collision energy above 13 GeV. The nonjet quadrupole 
centrality and pt dependence [13, Hl[ are independent of 
SS or AS jet structure. A unique characteristic of the 
nonjet quadrupole relative to the SS 2D peak is zero cur- 
vature on ?7a within the STAR TPC rj acceptance, with 
important implications for IS geometry and the possible 
structure of flows on 7/a- 



VI. PERIODIC PEAK ARRAYS ON AZIMUTH 

Because azimuth <p is a periodic variable any ID struc- 
ture on 4>A can be described exactly as a discrete Fourier 



7 



cosine series 

g(M = ^2 F m cos(m(f> A ). (8) 

m=0 

In principle, a factorized 2D structure f{r]A)g{4>A) could 
be so expressed by multiplying Eq. (|8) through by the 
j]a factor (see App. [C|l. However, representing some 2D 
structure or its ID projection by a few terms of a single 
ID Fourier series can be misleading. In this section we 
consider the Fourier series representation of a periodic 
Gaussian peak array on ID azimuth and projection of a 
SS 2D Gaussian onto azimuth. In the next section we dis- 
cuss possible confusion between jet-related Fourier com- 
ponents of the SS 2D jet peak and the nonjet quadrupole 
(see Apps. Eland IE). 

A. Fourier representation of ID peak arrays 

Because ID dihadron azimuth distributions are peri- 
odic, the peaks observed at <pA = (SS, same-side) and 
4>a = I" (AS, away-side) are actually elements of sepa- 
rate periodic peak arrays described by cosine series. The 
SS array is centered on even multiples of tt, the AS ar- 
ray on odd multiples. Nearest array elements outside a 
2-7T interval (image peaks) produce significant structure 
within the observed interval and should be included in 
fit models. 

Each peak array (SS or AS) can be represented by a 
Fourier series of the form 

oo 

S{<j>A',^(i> A ,n) = ^2 F m ^ n cos(m [(f) a ~ nir]), (9) 

m=0 

where the F m , n are functions of r.m.s. peak width oa^ 
defined below. Since n is even for SS peak arrays (+) and 
odd for AS arrays (— ) odd multipoles must be explicitly 
labeled as SS or AS. The terms represent 2m poles, e.g. 
dipole (to = 1), quadrupole (to = 2), sextupole (to = 3) 
and octupole (m = 4) , referring to cylindrical as opposed 
to spherical multipoles. 

Fig. [6] (left panel) illustrates the sum of peak arrays 
(solid points) for SS and AS peaks extending beyond one 
2ir period. The SS Gaussian peak array is the dash- 
dotted curve, the AS array with a A (/> ~ 7r/2 is the dashed 
curve (approximately pure dipole in this case) . The dot- 
ted curve is the quadrupole term of the SS peak array, 
which would add a large "nonflow" contribution to w|{2} 
inferred from that distribution. 

Fig. |S] (right panel) shows the Fourier amplitudes F m of 
a unit-amplitude Gaussian array for the first five terms 
(to e [0,4]) of Eq. jS) as functions of the r.m.s. peak 
width 

F m (a^) = V2A^Aexp(-TO 2 a^ A /2). (10) 

As peak width a c / >A increases, the number of signifi- 
cant terms in the series decreases. The limiting case is 




Gaussian width c^ a 

FIG. 6: Left: Periodic arrays of SS (dash-dotted) and AS 
(dashed) peaks. The SS peaks are Gaussians. The AS peaks 
are well-described by a dipole. The dotted sinusoid corre- 
sponds to the m = 2 Fourier component of the SS peaks. 
Right: Fourier amplitudes Fm of a unit-amplitude Gaussian 
[Eq. (1101) ] vs peak width . 

a 4>& ~ 7r /2, for which the peak array is approximated 
by a constant plus dipole term [1 + cos((/>a)]/2 (SS) or 
[1 — cos(</>a)]/2 (AS). For narrower (SS) peaks terms 
with to > 1 become significant, and a Gaussian func- 
tion is the more efficient representation. In particular, 
for (i0 A ss 0.65 (typical for the SS jet peak) jet-related 
quadrupole amplitude F2 ~ 0.22 represents the domi- 
nant jet-related nonflow contribution to v%{2} ~ v%{EP} 
data. 



B. Projecting a 2D Gaussian onto ID azimuth 

The relation between an SS 2D peak projected to a 
narrow SS ID Gaussian on azimuth and its Fourier com- 
ponents comes into play in Sees. I Villi and HXl which ad- 
dress recent claims of higher flow harmonics. The SS 2D 
peak, well-modeled by a 2D Gaussian, is distributed on 
(jia, </>a) within angular acceptance (A77, 2ir). We wish to 
determine the amplitude of the equivalent SS ID Gaus- 
sian projected onto azimuth c/)a- We first consider the 
case that the entire r\ acceptance A?y is projected onto 
4>A- We consider more-complex 77 exclusion cuts in Sec. IXl 
and App. |Dj The projection factor is given by 

GK - A ^ ) = CMAv-x) (1) 

= V^Cerf(l/C) - C 2 [1 - exp(-l/C 2 )] , 

with s = \/2cr I)A and C = s/Arj. G(a VA , Arf) —> 1 as £ — > 
00. We consider an application with SS peak parameters 
for 0-5% central Au-Au collisions and Ary = 2. 

Fig. [7] (left panel) shows Eq. (fTU) for ct^ a = 0.65, with 
the first few multipole coefficients marked for reference 
(open circles). Figure [7] (right panel) illustrates 2D^1D 
projection and Fourier decomposition of the SS peak. 
The SS (dash-dotted curve) and AS (dashed curve) peaks 
model jet correlations from 0-5% central Au-Au collisions 
where the nonjet quadrupole amplitude is consistent with 
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FIG. 7: Left: Evaluation of Eq. (fHJ)l for four values of m. 
Right: Typical results for ID azimuth correlations from near- 
central Au-Au collisions. Same-side peak (dash-dotted curve) 
and broad away-side peak approximated by a dipole (dashed 
sinusoid). The narrow SS peak can be decomposed into sev- 
eral sinusoids (m = 1...4, dotted curves). 

zero @, [20|. From Fig. [3] (v « 6) we obtain A 2D ~ 0.7, 
calculate G(2.25, 2) 0.85 from Eq. (JTTJ) and obtain the 
F m from Fig. [7] (left panel). The jet-related multipole 
amplitudes are then given by 

2p (b)v 2 m {SS}(b)^F m (^ A )G(a v& ,A v )A 2D (b).(12) 

Figure [7] (right panel) shows the corresponding azimuth 
multipoles as dotted sinusoids with amplitudes from 
Eq. Jill). Note that A Q {SS} = p (b)v${SS}(b) is 
the jet-related quadrupole amplitude, and 2Aq{SS} = 
0.225 x 0.85 x 0.7 = 0.135 defines the m = 2 sinusoid. 



VII. MINUETS vs NONJET QUADRUPOLE 

The interplay between (mini)jet structure and nonjet 
quadrupole contributions to measured v 2 data plays a 
central role in the interpretation of RHIC data. The non- 
jet quadrupole Aq{2D} is a unique phenomenon with 
centrality, energy and pt dependence distinct from mini- 
jets. Some v 2 analysis methods confuse jet and nonjet- 
quadrupole structure, leading to possible crosstalk be- 
tween jets and quadrupole in v 2 data (see App. [B]). We 
here focus on azimuth quadrupoles (m — 2) within a con- 
tiguous r\ acceptance (e.g., \f]\ < 1), then consider higher 
multipoles and r\ exclusion cuts in subsequent sections. 

A. Jet-related vs nonjet quadrupole 

The azimuth quadrupole measured by v 2 is conven- 
tionally attributed to "elliptic flow." Contributions to 
v 2 from possible nonhydro mechanisms are called "non- 
flow." In the present context we refer instead to a nonjet 
quadrupole (what might be attributed to elliptic flow) 
and a jet-related quadrupole (v 2 contribution mainly 
from jets and mainly from the SS 2D jet peak) . The dis- 
tinction between flow and nonflow has been extensively 
discussed (e.g., Ref. [35| and see App. [Bl. 



Figure [8] (left panel) shows the centrality dependence 
of SS 2D peak amplitude A 2 r> and three quadrupole am- 
plitudes related by A Q {2} = A Q \2D} + A Q {SS} [HJ. 
A Q {2D} is defined by Eq. © [H and An{SS} by 
Eq. (fT2j) using SS peak parameters from Ref. [4[ summa- 
rized in Fig. [31 Aq{2} (dotted curve) is then a prediction 
of v 2 measurements derived from the ID projection onto 
azimuth of all 2D angular correlation structure, corre- 
sponding to measured v 2 {2} ~ v 2 {EP} [3ol |40|. 




FIG. 8: Left: SS 2D (jet) peak amplitude A 2D , "nonflow" SS 
peak quadrupole component Aq{SS} and nonjet quadrupole 
A 2D amplitudes, with A Q {2} = Aq{2D} + A Q {SS} [H. 
Right: Quadrupole amplitudes Aq converted to conventional 
measure v 2 . Open squares are v 2 {2} data from [35l |. 

Figure [5] (right panel) shows ^{method} trends ob- 
tained from the corresponding Aq {method} curves in 
the left panel via Eq. ([5]). Also included are v 2 {2} data 
(open squares) from Ref. (35[. Good agreement between 
data and prediction (dotted curve) is evident. Thus, from 
2D nonjet quadrupole and minijet measurements we ac- 
curately predict v 2 {EP} ps v 2 {2} published data. The 
prediction does not include small contributions to i>2{2} 
from HBT and electron pairs (more significant for pe- 
ripheral A-A collisions) which are excluded from A 2 r> by 
the 2D model- fit procedure For statistically well- 
defined v 2 methods (e.g., v 2 {2} s» v 2 {EP}) the jet con- 
tribution can be estimated accurately and "flow" can be 
distinguish from "nonflow." This exercise for m = 2 is 
generalized to higher multipoles in following sections. 



B. ID ZYAM subtraction: m = 2 

Dihadron angular correlations are projections onto ID 
azimuth of 2D angular autocorrelations. Dihadron cor- 
relation analysis includes "trigger-associated" pt cuts in- 
tended to enhance jet structure relative to background. A 
nonjet combinatoric background is estimated by the zero 
yield at minimum or ZYAM procedure and subtracted 
from the sibling-pair density distribution to obtain an 
estimate of jet-related correlation structure. The back- 
ground estimate relies on published v 2 (p t ) data [26l[35j. 

Two problems arise from ZYAM subtraction: (a) the 
ZYAM offset estimate is not valid for overlapping peaks 
(as encountered in more-central A-A collisions) leading 
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to large errors in the apparent zero offset and inferred 
peak shapes and (b) published U2(pt){ me th°d} data may 
include substantial contributions from the SS 2D jet peak 
(nonflow). The result is underestimation of jet fragment 
yields and distortion of inferred jet correlations [27[. 
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FIG. 9: Left: Simulated dihadron correlation data (bold solid 
curve) for central 6 = Au-Au collisions relative to true base- 
line, with SS ID Gaussian (dash-dotted curve) and AS dipole 
(dashed curve). Dotted curves labeled D, Q and S are re- 
spectively the dipole, quadrupole and sextupole Fourier com- 
ponents of the SS ID Gaussian. Right: Original correlation 
data (bold solid curve, total jets) and bold dashed curve il- 
lustrating the result of ZYAM subtraction with biased back- 
ground V2 including 50% of the SS ID Gaussian (jet-related) 
quadrupole component. The dash-dotted curve illustrates the 
origin of the AS structure (see text). 



dash-dotted curve represents the sum S/2 + net dipole, 
demonstrating the origin of ZYAM-induced distortion of 
AS structure. Such distorted and suppressed jet man- 
ifestations are used to support claims for formation of 
an opaque bulk medium with anomalous properties in 
RHIC collisions. The distortions arising from v 2 over- 
subtraction in ZYAM analysis are worsened by newly- 
introduced analysis methods described in the following 
two sections. 



VIII. TRIANGULAR FLOW 

In a follow-up to ZYAM subtraction the SS sextupole 
component "released" from the SS ID Gaussian by v 2 
oversubtraction is redefined as the FS manifestation of a 
"triangularity" component of IS geometry coupled to ra- 
dial flow Q. The m = 3 Fourier component of the SS 2D 
peak is then identified as "triangular flow." The strat- 
egy eliminates remaining manifestations of the SS 2D jet 
peak as such in the ID projection, in part by suppressing 
the critical tja curvatures which distinguish SS jet struc- 
ture from nonjet structure. In this and the following 
section we summarize recent arguments favoring "higher 
harmonic flows" V m in the final state. In Sec. [X] we re- 
spond with V m predictions based on a two-component 
minijet-quadrupole model of FS correlation structure. 



Figure |pj (left panel) shows simulated data from 200 
GeV central (b = 0) Au-Au collisions (from Sec. [V} 
projected onto azimuth (bold solid curve). For central 
Au-Au collisions the nonjet quadrupole is consistent with 
zero [2(| ■ Correlation structure consists entirely of the SS 
ID Gaussian (dash-dotted curve) and AS dipole (dashed 
curve). The zero offset is well-defined based on model 
fits to 2D histograms. Also shown are SS ID peak (jet- 
related) multipoles D, Q, S (dipole, quadrupole, sex- 
tupole). In 2D angular correlations the SS multipoles 
with their large curvatures on t/a are distinct from the 
AS dipole and nonjet quadrupole which have negligible 
curvatures within \rj\ < 1 (see App. ICj) . 

Fig. [5J (right panel) shows the result of ZYAM back- 
ground subtraction (dashed curve). Half the SS ID peak 
quadrupole amplitude V2{SS} is typically included in the 
v-i used for ZYAM subtraction (nonflow bias) 27]. The 
SS ID Gaussian is then effectively divided into a half- 
amplitude Gaussian and "liberated" SS dipole and sex- 
tupole terms. The original data (total jets, bold solid 
curve) and a half-amplitude SS ID peak (SS/2, lower 
solid curve) are shown for comparison. What remains 
after subtraction (dashed curve) is the sum of the half- 
amplitude SS Gaussian (SS/2), a half-amplitude sex- 
tupole (S/2) and a small-amplitude (negative) net dipole 
(from near cancelation of SS/2 + AS dipoles). The small 
net (negative) dipole is reinterpreted as "bulk momen- 
tum conservation," and the S/2 SS sextupole produces 
relatively large distortions in the surviving AS struc- 
ture interpreted as evidence for "Mach cones" [4l| . The 



A. Example: ID Fourier analysis 

Figure [10] is similar to the right panel of Fig. 1 from 
Ref. [9( which includes an analysis of STAR dihadron 
data. The points in Fig. QJJ are simulated 10-20% cen- 
tral 200 GeV Au-Au data derived from a STAR data 
analysis jij. The corresponding 2D histogram is similar 
to that shown in Fig. [5] (c) with a somewhat different 
centrality (10-20% instead of 20-30%). The 10-20% pa- 
rameters from Eq. (0} arc A 2 d = 0.77, 2A Q = 0.18 and 
Ad = 0.58. The amplitude of the SS ID Gaussian pro- 
jected from 2D is A 1D = 0.85 x 0.77 = 0.65. Projection 
factor G = 0.85 (with rj exclusion cut, see App. [D]) cor- 
responds approximately to the restricted ( "long range" ) 
r] interval employed in Ref. [§] {-q cut noted in Fig. [TU]) . 

The data points shown in Ref. @ represent pair ra- 
tio (r) = p/p re f and are plotted on one quadrant of 
(^Aj^a)- Thus, conversion of parameters Aio, Aq and 
Arj is required. The conversion factor from Eq. (Q} to 
Ap/p ref = p/p^f - 1 is 1/poO), with Po(10-20%) w 70. 
An additional factor 4 is required to convert from a single 
quadrant to the full angular acceptance. For the STAR 
data used in Ref. [t| a pt > 0.8 GeV/c cut was applied, 
which results in a somewhat narrower (on azimuth) SS 
2D peak with reduced amplitude. 

Figure [lOj also shows model elements for the SS ID 
Gaussian and AS dipole (dashed curves) and their sum 
(dash-dotted curve labeled "jets") with the true baseline 
established by the 2D fit to data. Fourier components of 
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"long-range" (|^a| > 1) and "short-range" (|?7a| < 1) 
projections is minor. Short-range correlations may re- 
veal somewhat larger "triangular flow." There is no ac- 
knowledgment of possible jet structure in minimum-bias 
angular correlations. 

In analysis based on ID projections the large curvature 
on r)A of the SS 2D peak is suppressed. But that curva- 
ture distinguishes the SS 2D peak from the AS dipole 
and nonjet quadrupole, both approximately uniform on 
7/A within the relevant 7/ acceptance. Thus, the sums 
2AQ+F2A10 and 2FiAi£,—Ao combine numbers belong- 
ing to different categories of 2D structure with misleading 
results. The Fourier components of the projected SS 2D 
peak cannot be considered individually (e.g., "triangular 
flow" ) because they must sum to zero in the AS azimuth 
region to describe measured 2D data. This approach 
marks a reversal of conventional attempts to discrimi- 
nate "elliptic flow" (nonjet quadruple) from "nonflow" 
(jet- related quadrupole), for example with four-particle 
«2{4} or Lee- Yang zeroes methods [H[ (see App. IB")) . 



FIG. 10: Fourier analysis of azimuth correlations for 10-20% 
central 200 GeV Au-Au collisions. Simulated data based on 
Ref. 0] are shown as points projected from a 2D histogram. 
Jet structure inferred from 2D fits is shown as dashed and 
dash-dotted curves. Solid curves represent Fourier compo- 
nents similar to those derived in Ref. |J but inferred in this 
case from fitted jet structure and nonjet quadrupole (see text). 
The curve labeled m = 4 can be compared with the fit resid- 
uals in Ref. [§]. 



the ID data were derived in Ref. [9j as follows. A four- 
term Fourier series (m = 0-3) was fitted to the data to 
obtain amplitudes for dipole, quadrupole and sextupolc 
terms. The Fourier fit was subtracted from the data to 
obtain residuals dominated by an m = 4 octupole com- 
ponent. 

Results from such a procedure can be anticipated from 
measured properties of the SS 2D Gaussian, AS dipole 
and nonjet quadrupole as follows. From Fourier coeffi- 
cients F m (cr0 A ), with <70 A = 0.65 for the SS ID peak as 
determined in Sec. I VII we obtain total quadrupole ampli- 
tude 4(2A Q +F 2 A 1D )/70 = 0.019 ("flow" + "nonflow" or 
nonjet + SS quadrupoles respectively), total dipole am- 
plitude A{2F 1 A 1D - A D )/7Q = -0.0017 (SS + AS dipoles 
respectively) and SS sextupole 4F 3 A 1D /70 = 0.0029. 
Residuals should be dominated by SS octupole amplitude 
AF4A1D/7O = 0.0007. The corresponding solid curves in 
Fig.[TU]are consistent with the results presented in Fig. 1 
of Ref. @. 

Reference @ interprets results from the Fourier analy- 
sis (solid curves) to indicate that there is a large "elliptic 
flow" component (m = 2), a significant "triangular flow" 
component (m = 3) and a small negative dipole (m = 1) 
representing "global momentum conservation." The de- 
scription refers to "long-range" correlations, but exami- 
nation of Fig. [2] (c) reveals that the difference between 



B. Triangular flow centrality systematics 

Figure [TT] (left panel) shows ^{method} plotted in 
the format of Ref. The curves are based on STAR 
minimum-bias angular correlations presented in Sec. [V] 
The SS 2D peak 77 width plays the dominant role in 
systematic variations. This plot is a variant of Fig. [5] 
(right panel), but the 77 exclusion cut applied in this case 
(I^aI £ [2,4]) insures much stronger "extinction" of the 
SS 2D jet peak contribution in more-peripheral A-A colli- 
sions (below the sharp transition at v s» 3 or N part < 50) 
than the cut |?7a| € [1, 2] applied in Fig. [Poland Sec. IX CI 
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FIG. 11: Left: Comparison of three V2 methods obtained from 
the parametrizations in Sec. [V] with rj exclusion cut |?7a|£ 
[2, 4] and plotted on N part as in Fig. 5 (left panel) of Ref. [Sj]. 
Right: Ratio of sextupole to quadruple amplitudes vs N par t 
comparable to Fig. 8 (left panel) of Ref. [I| 

Figure [TTJ (right panel) shows ratio v${SS}/v%{2} 
which includes contributions from the SS peak in both 
numerator and denominator but contributions from non- 
jet quadrupole t>|{2.D} in the denominator only. The 
ratio also falls to zero below N part = 50, because the 
only contribution to numerator v^{SS} is from the SS 
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2D peak. With increasing centrality tj|{2Z?} in the de- 
nominator falls to zero and the ratio must increase to the 
limit 3F3/F2 « 1 consistent with 2D correlation data. 
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FIG. 12: Left: Sextupole amplitude obtained from the 
parametrizations in Sec. [V] with n exclusion cut \t)a\ £ Pi4] 
and plotted on N part as in Fig. 5 (right panel) of Ref. [jj]. 
Right: The same result plotted on participant path length v, 
showing the correspondence to the sharp transition in SS 2D 
peak properties in Fig. [3j especially the n elongation. 

Figure [T2l shows vl{SS} for two plotting' formats com- 
parable to Fig. 5 (right panel) of Ref. Q. Note that 



,{2} = v 2 m {SS} for 



3. The sextupole compo- 



nent becomes nonzero only above the sharp transition 
in SS 2D peak properties where the SS peak elongates 
on i] into the "ridge-like" tja acceptance. Per-pair mea- 
sure v%{SS} decreases strongly above the sharp transi- 
tion at v sa 3. In contrast, per-particle sextupole ampli- 
tude Asl^S} would continue to increase monotonically 
with centrality, reflecting increased jet production. 

The results from this demonstration based on 
previously-measured minijet 2D angular correlations are 
generally duplicated by results from Ref. @. Reported 
"triangular flow" measurements are consistent with the 
properties of a monolithic SS 2D jet peak. Further de- 
tails are presented in Sec. |Xj The effects of 77 exclusion 
cuts are considered in detail in Sec. IX CI 



IX. LONG-RANGE CORRELATIONS 

In a followupto triangular flow Ref. [lOj disputes a 
claim in Ref. |42J that dependence on the reaction-plane 
angle of ZYAM-subtracted long-range (on 77) azimuth 
correlations signals nonflow. Instead, all such "ridge- 
like" correlations should be attributed to collective flows, 
and ZYAM-subtracted data are said to be consistent with 
that hypothesis: "...the measured dihadron correlation at 
large A77 consists almost entirely of the lowest few Fourier 
components, each of which can be quantitatively under- 
stood as coming from collective flow (plus global momen- 
tum conservation)..." (Io| . Long-range azimuth correla- 
tions are said to originate at early times and therefore 
must represent the "collective behavior of the system" 
(and global momentum conservation). Here we present 
evidence contradicting claims in Ref. [Hj] and Ref. (Toj . 



A. Fourier series analysis of long-range structure 

A Fourier series analysis of unsubtracted "ridge-like" 
trigger-associated ID azimuth correlations (with trig- 
ger related to the event plane) is shown in Fig. 1 of 
Ref. [lOj]. 77 exclusion cuts 0.7 < \t]a\ < 2 reject a 
short-range or jet-like part of the pair acceptance on 
?7a in favor of a long-range or ridge-like part of the ac- 
ceptance. The analysis reports four significant Fourier 
coefficients Vm{2} (m = 1-4), where V m = ujj,^ is 
nominally a product of trigger and associated ua(pt) val- 
ues. The V<i{2} quadrupole term dominates mid-central 
Au-Au collisions. Absence of higher harmonics (to > 4) is 
considered remarkable but should be expected given the 
widths of peaked structures observed by previous analy- 
ses [1, 13,113; as discussed in Sec. IVII 

Inferred quadrupole and octupole terms V^{2} and 
V^{2} are attributed solely to elliptic flow. The am- 
plitudes are said to vary relative to event-plane angle 
as cos(20 s ) and cos(40 s ) as expected for flows, with 
<f> s = \(f>t — iPep\ (trigger and event-plane angles). Fitted 
Va{2} values are said to exceed "measured" V2{ZA} data 
(estimated V2 [42]) by a significant amount (attributed 
to nonflow). But since the nonflow appears to have the 
same dependence on p t and <fi s , Vi{ZA\ must underes- 
timate the true flow, and to = 2,4 "nonflow" must be 
flow. 

Conjectured flow terms Vl{2} and ^{2} are not ex- 
pected to correlate with the reaction plane, and the 
Fourier analysis in Ref. [l(| is interpreted to confirm 
that expectation. The to = 3 sextupole term increases 
monotonically with p t . It is therefore concluded that the 
to = 1,3 terms are also consistent with flows. Thus, all 
long-range or "ridge-like" structure must be flows. 

However, the "nonflow" components of V^{2} and 
I4{2} (excesses over corresponding ZA estimates) do not 
actually follow the sinusoid trends expected for flows (see 
following subsections). V3 is said be independent of <p s , 
implying a flow interpretation, but the V3 data from four 
of five pt intervals show significant increases, and real 
"nonflow" (jets) may or may not have a significant de- 
pendence on cj) s . There is insufficient sensitivity to details 
of pt dependence that might distinguish "flows" from jet 
structure. Thus, the results in Fig. 1 of Ref. [l(|, ex- 
amined sufficiently differentially, are either inconclusive 
or actually confirm the presence of nonflow (jets) in the 
dihadron data. In the rest of this section we examine the 
case for m = 2 in detail. 



B. Event-plane-related ZYAM subtraction 

The ZYAM discussion in Ref. [13] is based on an anal- 
ysis of dihadron correlations which included free fits with 
a model function applied to dihadron data as a check on 
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the ZYAM subtraction 0. The fit model 



dN v 



= B {1 + 2V 2 cos(2</> A ) + 2T/ 4 cos(4^ A )}(13) 
+ Aice-^OW^J 2 + Abcos(0a -tt) 



is similar to Eq. Q of the present analysis. Based on 
the fit results Ref. 42j concluded that free model fits 
must be incorrect for two reasons: (a) "In order to elimi- 
nate the away-side double-peak... flow [fitted V2, V4] that 
is much larger than experimentally measured [V m {ZA} 
defined below] is required...." and (b) "...deviations of 
the fitted flow modulations from the measured ones vary 
slice to slice [cj> s bins], which should not be the case if 
the measured flow parameters... were simply in error" [em- 
phasis added]. Thus, "...the fit model [Eq. (fT5])]...cannot 
be the correct functional form to describe the dihadron 
correlation signal...." But that conclusion doesn't fol- 
low from the fitting exercise, which actually demonstrates 
that the unsubtracted dihadron data are well described 
by Eq. (H3J). 

The values of V2{ZA} used to define the background 
for event-plane-related ZYAM subtraction are derived 
from Eq. 4 of Ref. [12], which assumes prior knowledge 

of the product (ptjv^ (pt) for trigger and associated 
Pt bins and the event-plane resolution from separate nu- 
merical analysis. There is reason to question such data 
based on independent 2D correlation analysis 21 1. 

Figure [T3] is similar to Fig. 1 (left column, second 
panel) of Ref. [HJ. The figure shows "measured" V 2 {ZA} 
valu es ( ZYAM analysis — > ZA) reported in Table II of 
Ref. [42| (open squares, note error bars within) based on 
parameters in Table I. The dash-dotted curve is Eq. 4 of 
Ref. [42[ with the same parameters from Table I inserted, 
confirming consistency. The "fitted" V 2 {1D} values in- 
ferred from ID free fits (reconstructed here from correc- 
tion factors provided in Ref. (42[) are shown as open cir- 
cles. The free-fit results also constitute measurements, 
but derived directly from fits to the dihadron data in 
question rather than indirectly from separate analysis of 
other particle data. 

Variation of correction factors "from slice to slice" (a) 
may be irregular due to the rapidly changing magnitude 
(and sign) of V%, but the vertical offset from one data 
trend to another varies slowly and smoothly. "Much 
larger" V2 values (b) are not observed. Absolute changes 
in inferred quadrupole amplitudes are actually modest, 
within typical systematic uncertainties for v 2 {EP\ data 
(see below). The solid curve is Eq. (4) of Ref. [42j with 
the event-plane resolution increased by 10% and the trig- 
ger v 2 reduced from 0.16 to 0.1, consistent with Ref. |2l| . 

It is important to note that the small difference be- 
tween the solid curve and the dash-dotted curve in Fig. [1^1 
accounts entirely for the difference between jet structure 
undistorted and increasing in amplitude from in-plane to 
out-of-plane and jet structure increasingly distorted and 
suppressed from in-plane to out-of-plane from Ref. [42j . 
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FIG. 13: Measured quadrupole amplitudes from two methods. 
V^-fl-D} data (open circles) are derived from free model fits 
to ID azimuth correlations, and V2{ZA} data (open squares) 
are estimated from published V2 data [4^]. Error bars are 
shown within the open symbols. The V2{2} trend (dashed 
curve) is based on an estimate of the SS 2D peak contribution 
V2{SS} summed in quadrature with V2{1D}. Thin vertical 
lines mark zeros of cos(20 s ) and cos(40 s ). The dash-dotted 
and solid curves are both obtained from Eq. (4) of Ref. [42l ] 
(see text). 



C. Argument for reinterpretation of ZYAM results 

ZYAM subtraction is intended to remove "elliptic flow" 
and an uncorrelated combinatoric back grou nd from di- 
hadron azimuth correlations. Reference [42| argues that 
since ZYAM-subtracted "ridge-like" correlations depend 
on reaction-plane angle <p s (e.g., increasing suppression of 
jet-like structure) the surviving structure must represent 
"nonflow," since recently-proposed higher odd flow har- 
monics (e.g., V3) should not relate to the reaction plane. 

Based on its own inferred <p s and pt trends Ref. [l(| 
responds that additional flow structure must remain af- 
ter ZYAM subtraction, the amount depending on <fi s . 
ZYAM v m "should come from an independent measure- 
ment that does not contain a contribution from non-flow 
correlations.... Such a measurement does not exist." Be- 
cause the v m are not known "...the result [of ZYAM sub- 
traction] will have significant contributions from flow." 
In other words, unknown systematic v m biases must be 
such that measured v m always underestimate "real" flows 
by overestimating nonflow. The upper limit estimated 
by V 2 {ZA} is therefore too low. One should impose a 
Fourier analysis on all long-range or "ridge-like" correla- 
tions to obtain F m {2} representing "real" flows. With 
that definition no nonflow (jets) could survive "back- 
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ground" subtraction — all jet structure must be redefined 
as flows. 



D. Estimating the nonflow contribution 

In Ref. [l(| we identify the term "estimated flow con- 
tribution" with V2{ZA} and "extracted Fourier compo- 
nent" with fitted T^{2} ~ V 2 {EP}, which is said to mea- 
sure "real" flow. The difference V 2 {2} - V 2 {ZA} is de- 
fined as "nonflow," which is claimed to be overestimated 
based on inferred 4> s and p t trends. Reference [l(| con- 
cludes that nonflow should be reduced to zero and V2{2} 
adopted as "real" flow to be used in ZYAM subtraction, 
along with other higher harmonics similarly determined. 
The concept is illustrated in Fig. 2 of Ref. [10], with con- 
trasting subtraction results from V 2 {ZA} (left panels) 
and from V2{2} (right panels). 

Based on study of 2D correlations and their ID projec- 
tions we find that V 2 {2} = V 2 {SS} + V 2 {1D} ^M,m- 
And V 2 {SS} is quantitatively predicted by the SS 2D 
peak with its large curvature on 7/a, whereas V2{1Z?} ~ 
V 2 {2D} represent an azimuth sinusoid with no curvature 
on 77 a- We can therefore predict V2{2} from Ref. [l(| 
based on the free-fit results and the quadrupole compo- 
nent of the measured SS 2D jet peak (nonflow) as fol- 
lows. The SS peaks derived from free fits and shown 
in Fig. 4 of Ref. [42j (dashed curves) increase in ampli- 
tude with (f> s smoothly over the interval 0.7-1.2. The 
background constant for the 1-2 GeV/c p t cut interval is 
B w 3.8. Given Fourier coefficient F 2 0.2 we obtain 
for the SS peak quadrupole 2V 2 {SS} = 0.2(0.7-1.2)/3.8 
or V 2 {SS} 0.02-0.03. Combined with V 2 {1D} from 
the free fit (solid curve) in Fig. [13] we obtain V2{2} 
as the dashed curve, which is consistent with Fig. 1 of 
Ref. [To| | ■ Positive-definite nonflow V 2 {SS} (jet-related 
quadrupole) increases monotonically with <fi s , a trend in- 
consistent with flow expectations. 

Given the structure of 2D angular correlations, non- 
flow relevant to ZYAM is more properly defined as 
V 2 {ZA} - V 2 {1D}, the excess of ZYAM v 2 over the non- 
jet (fitted) quadrupole value. The nonflow upper limit is 
then V 2 {2} - V 2 {ID}. Subtracting V 2 {2} would achieve 
complete removal of the SS peak from projected ID di- 
hadron correlations. But the SS peak remains the most 
prominent structure in 2D angular correlations and can- 
not be removed by such ID subtraction procedures. Re- 
duction of nonflow bias to zero in the V 2 estimate is ac- 
tually achieved by invoking V2{1Z?} from the free fits. 



X. PREDICTING v m MEASUREMENTS 



Nongraphical numerical methods (e.g., conventional v m 
measurements) and dihadron correlations depend explic- 
itly or implicitly on ID projection followed by a single 
Fourier-series decomposition of the azimuth projection. 

We generalize the "flow vs nonflow" case for m = 2 
to higher m by projecting all measured 2D structure to 
ID and invoking a single Fourier series to predict recent 
v m {2} measurements. We introduce 77 exclusion cuts 
to simulate attempts to distinguish (remove) jet struc- 
ture (nonflow) from conjectured hydro phenomena. This 
treatment emphasizes ID projections onto azimuth. Ex- 
tension to 2D correlations is described in App. [Cj 



A. ?7a dependence of inferred v m 

Measurements of the 77 dependence of "triangular flow" 
and other higher multipoles can be simply predicted from 
the 2D correlation model of Sec. [V] Since all higher mul- 
tipoles are derived solely from the SS 2D peak, its pa- 
rameters determine all multipole systematics. 




FIG. 14: (Color online) Left: Simulated angular-correlation 
data from 0-5% central 200 GeV Au-Au collisions. An AS 
dipole component has been removed. The remainder is the 
SS 2D peak well-described by a 2D Gaussian. Right: The 
sextupole component of the SS 2D Gaussian from the left 
panel determined as a function of 77A , predicting what would 
result from such a "triangular flow" V3 analysis applied to the 
same particle data. 

Figure [14] (left panel) shows a correlation model of the 
SS 2D peak for 0-5% central 200 GeV Au-Au collisions. 
It is the histogram in Fig. [2] (d) with the AS dipole re- 
moved. Fig. 03] (right panel) shows the corresponding 
vj{SS}(n A ) defined by 



2p (b)vl{SS}(r] A ) = F 3 (a 0A )A 2D exp 



Va 
2< 



(14) 



Note that the mean value of that trend is consistent with 
the value of for most-central collisions from ID 

projection onto azimuth in Fig. 1151 (right panel). 



Based on model fits to 2D histograms, including mini- 
jet and nonjet-quadrupole structure, we can predict any 
v m measurement for any 77 exclusion cuts. In contrast, ID 
projections onto 77A or cf> A cannot predict 2D structure 
because the projections abandon essential information. 



B. v m centrality trends in contiguous 77 acceptance 

The centrality dependence of any multipole integrated 
over a contiguous 77 acceptance symmetric about the ori- 
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gin denoted by A77 can be predicted from the 2D data 
model described in Sec. [V] The multipoles derived from 
the SS 2D peak are given by Eq. (jl2]l . 




80 100 



% Central % Central 

FIG. 15: Left: Quadrupole components from 200 GeV Au-Au 
data vs centrality satisfying the relation t>2{2} = V2{2D} + 
v%{SS} [2J], where v%{SS} is the jet-related quadrupole com- 
ponent of the SS 2D peak. Right: Corresponding higher mul- 
tipoles of the SS 2D peak, predicting the result of "higher 
harmonic flow" analysis applied to the same particle data. 

Figure [15] shows quadrupole, sextupole and octupole 
amplitude predictions in the form ^{method} (b) for 200 
GeV Au-Au collisions, where SS (dash-dotted curves) in- 
dicates a multipole derived from the SS 2D peak, and 2D 
(solid curve) refers to the nonjet quadrupole inferred from 
2D fits to data (Sec. IV Cj) . All SS multipoles are predicted 
to have the same rj& and centrality dependence, which 
would then confirm a common source (SS 2D peak). The 
common SS trends are very different from the unique 
nonjet quadrupole u|{2Z3} contribution. The good agree- 
ment between published data (open squares) and the 
predicted trend (dashed curve) is notable, since 

the prediction combines two independent measurement 
programs ( minij et systematics [H and nonjet quadrupole 
systematics |20() unrelated to the published data analy- 
sis m. 

C. v m centrality trends with n exclusion cuts 

It is argued that SS angular correlations can be sep- 
arated by cuts on r\ into a "short-range" or "jet-like" 
region at smaller 77A and a "long-range" or "ridge-like" 
region at larger 77 a [9]. Structure in the ridge-like region 
is attributed to nonjet mechanisms, for instance initial- 
state geometry coupled with radial flow. We derive an 
expression for the effect of r\ exclusion cuts in App. [D] 
assuming that particle pairs are formed from symmetric 
disjoint r\ intervals [ — 772 , — ?7i] and [771 , 772] to obtain pro- 
jection factor G m (a VA ; 771, 772). Multipole amplitudes are 
then defined in terms of 2D model elements by 

2 Po(b) v'£{SS} = F m (<70 A )G m (o-^ A ;7/i,7/2)A 2 D(b), (15) 

the prime indicating an 77 exclusion cut. 

Figure [TO] shows multipoles from the full 77 acceptance, 
as in the previous subsection (upper curve for each line 
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FIG. 16: Illustration of "nonflow" reduction by 77 cuts. The 
bold solid curve is the parametrized nonjet quadrupole trend 
from Ref. [201 ]. For each pair of remaining curves the upper 
curve corresponds to the full pair acceptance \va\ < 2 (as in 
Fig. I15|) whereas the lower curve of each pair corresponds to 
a reduced pair acceptance 1 < \t)a\ < 2 intended to reduce 
nonflow. The v m {SS} trends are determined by the SS 2D 
peak properties presented in Fig. [3] with and without rj cuts. 



style), and from 77 exclusion cuts with 771 = 0.5 and 772 = 1 
(lower curve for each line style) which exclude the "short- 
range" interval |?7a| < 1. The solid curve shows nonjet 
quadrupole V2{2D} which is independent of such cuts. 
The trends reflect smooth elongation of the SS 2D peak 
on ?7a with increasing centrality, modulo a sharp transi- 
tion in SS peak properties just below v = 3. 

The p t -integral SS 2D peak shows no indication of be- 
ing composite (e.g., distinct jet and ridge components on 
?7a) for any Au-Au centrality. The SS peak is described 
in all cases by a single 2D Gaussian. The difference be- 
tween v m and v' m is entirely due to the changing SS peak 
?7A width relative to the special 77 cut. More of the SS 
peak is accepted by the 77 exclusion cut in more-central 
collisions. The apparent extent of SS peak rejection (ex- 
tinction) for more-peripheral collisions is minimized by 
plotting square root v' m rather than v!^ which measures 
the actual correlated pair number. Aside from 77 elon- 
gation all aspects of the SS 2D peak conform to pQCD 
jet expectations, including modification of fragmentation 
functions |7|. Note that in Figs. [TT1 and IT21 the 77 exclu- 
sion cut is extended to 2 < |?7a| < 4 leading to more 
complete extinction of the SS peak contribution for more- 
peripheral A-A collisions. Response of the V m {SS} to 77 
exclusion cuts is a very indirect way to study the evolving 
77 structure of the SS 2D peak. 
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D. Comparison with recent LHC results 

A test of v m predictions from the 200 GeV 2D data 
model is provided by comparison with a recent LHC anal- 
ysis which infers "higher harmonic flows" from Pb-Pb 
data Figure [T7l shows v m {2} measurements from 

2.76 TeV Pb-Pb collisions (points from Fig. 1 of Ref. [3) 
compared to the 200 GeV results summarized in Fig. \W\ 
(curves). The data from Ref. [l]| represent coefficients 
from a Fourier-series fit to the sum of all two-particle 
angular correlations projected onto ID azimuth with 77 
exclusion cut 1 < 1 77^. | < 2, denoted i> m {2,r; cut}. 
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FIG. 17: (Color online) The data points are from Ref. [13]. 
The various curves are obtained from Fig. [16] with <j^ a 
reduced from 0.65 to 0.60 and with an overall multiplier 
1.3. The relative magnitudes and centrality variations agree 
closely between data and curves, suggesting that the SS 2D 
peak at 2.76 TeV has similar properties to that at 200 GeV, 
as also noted in Ref. [43| . 

The curves in Fig. [17] are derived from Fig. [16] as fol- 
lows. The SS peak azimuth width C0 A is reduced from 
0.65 (200 GeV) to 0.60 (2.76 TeV). Azimuth width reduc- 
tion is consistent with the energy trend from 62 to 200 
GeV [4| , increasing the V3 and V4 trends slightly relative 
to v%. Given the SS peak width adjustment all curves in 
Fig- US] are multiplied by a common factor 1.3, the same 
factor attributed in Ref. [44J to an increase in the inclu- 
sive spectrum mean p t at the higher energy. For each 
line style the upper curve corresponds to a contiguous 
|?7a| < 2 acceptance while the lower curve corresponds to 
the nonflow exclusion cut 1 < |7?a| < 2. 

Given those adjustments the predictions derived from 
2D angular correlations at 200 GeV describe the Pb-Pb 
data at 2.76 TeV very well. The v 2 {2} trend follows the 
expectation for applied 7/ exclusion cuts, suggesting per- 



sistent presence of a sharp transition [4] in the SS 2D 
peak 77A width at the higher energy. The lack of V3 and 
t>4 data for more-peripheral centralities is unfortunate, 
because such data might provide a more direct test of SS 
2D peak systematics at the higher energy. The LHC data 
are also consistent with the t>2{2£>} nonjet quadrupole 
trend (solid curve) which shows no sensitivity to r\ exclu- 
sion cuts because the corresponding correlation structure 
is uniform on t/a within the 77 acceptance. 



XI. DISCUSSION 

This article addresses a prominent issue in the study of 
angular correlations at the RHIC: do jets or flows domi- 
nate the dynamics of high-energy heavy ion collisions? A 
trend has emerged recently to reinterpret jet angular cor- 
relations by projecting some part of the 77 acceptance onto 
ID azimuth and expressing the projection as a Fourier se- 
ries. The series elements may then be interpreted as flows 
(or global momentum conservation). The FS multipoles 
are in turn attributed to conjectured IS geometry fluctu- 
ations coupled to radial flow. We are obliged to question 
the validity of isolating individual Fourier components of 
a single 2D peak structure, no matter what its interpre- 
tation. We should also question hydro interpretation of 
any multipole, given measured systematics of the nonjet 
quadrupole, AS dipole and SS 2D peak [2fjj . 



A. Jet-related 2D angular correlation structure 

Correlations attributable to minimum-bias jets (mini- 
jets) appear as both angular correlations on (77A, 4>a) and 
momentum correlations on (p t ,pt)- The pQCD expec- 
tation for jet angular correlations is a nominally sym- 
metric same-side peak centered at the angular origin 
(intra-jet correlations) and a broad (on azimuth) away- 
side ID ridge (inter-jet correlations from back-to-back 
jets) peaked at </>a = fl". Jet correlations on {pt,Pt) 
(hard component) from p-p collisions are observed to 
be peaked near p t = 1 GeV/c and well resolved from 
a soft com pon ent which does not extend above 0.5 
GeV/c 0, 1H, HI, HJ. In P~P and peripheral Au-Au 
collisions the minijet hypothesis describes all data well. 
In more-central Au-Au collisions jet-like angular corre- 
lations are modified 0,01; but the jet hypothesis is cer- 
tainly not excluded by the data. 

In the context of hydro models the elongated SS 2D 
peak ("soft ridge") has been attributed to radial flow 
coupled to conjectured elongated structure in the initial- 
state transverse configuration space. However, there is 
presently no evidence to compel preference of flow con- 
jectures over a jet mechanism and substantial evidence 
that disfavors them [4^, [46| . 
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B. Jet-related terminology and partitions 

Misleading terminology may produce confusion regard- 
ing jets. Problems arise from arbitrary partition of kine- 
matic variables pt and rj based on conjectures about jet 
properties. Reference is made to "high-pt jets" or "semi- 
hard fragments." Hadrons are described as "soft" or 
"hard" based on partition of the pt axis. Low-p t ("soft") 
hadrons (e.g., below 2 GeV/c) are excluded from jets, 
reserved instead for "bulk" and "flow" phenomena. But 
for jets of any energy the majority of hadron fragments 
fall below 2 GeV/c 0, E3]. The terms soft and hard 
should describe IS momentum transfers between projec- 
tiles (hadrons or partons), not FS hadron fragments. 

The pair rj acceptance is arbitrarily divided into "short- 
range" and "long-range" regions, the former reserved for 
jet interpretations, the latter associated with conjectured 
bulk collective phenomena based on questionable causal- 
ity arguments. Long-range correlations are attributed 
to "early-time" interactions and are then said to be a 
manifestation of collective motion best described by the 
lowest few azimuth multipoles or "harmonic flows." But 
r\ is a measure of polar angle, not momentum per se. Jet 
structure at large r\ is not precluded. 

Biased terminology and unjustified partitions may in- 
correctly consign jets to a small fraction of the final state, 
nominally within |t?a| < 1 (or less) and p t > 2 GeV/c (or 
greater). The SS peak in more-central A- A collisions is 
arbitrarily partitioned into jet-like and ridge-like compo- 
nents with distinct production mechanisms. Any devia- 
tion of jet morphology from a notional ideal is attributed 
to nonjet (flow) mechanisms. Projection of some fraction 
or all of the SS peak structure onto ID azimuth leads to 
(mis) attribution of jet correlations to flows via Fourier 
analysis. 



C. Minimum-bias vs pt-conditional jets 

Factorization of two-particle momentum space into 
2D subspaces {pt,Pt) or (yt,yt) (transverse momentum 
or rapidity correlations) and (t^a^a) (angular correla- 
tions) retains almost all correlation information 0, HH ■ 
Minimum-bias angular correlations correspond to an in- 
tegral over the entire (pt,Pt) space. So-called "trigger- 
associated" correlations correspond to cut conditions 
(typically asymmetric) imposed on (j> t ,pt)- 

Isolation of nominal jet structure from nonjet struc- 
ture and combinatorial background typically follows dif- 
ferent strategies in the two cases. In minimum-bias anal- 
ysis a sibling-pair distribution containing correlations of 
interest is compared to a mixed-pair reference which is 
nominally uncorrelated. The net correlation structure is 
analyzed by 2D model fits which establish a zero offset 
from the fit, as described in Sec.fV] For pt-conditional or 
trigger-associated analysis the sibling-pair distribution is 
typically projected onto ID azimuth difference. A ZYAM 
background, estimated by the ZYAM offset criterion and 



published vi data, is subtracted from the sibling-pair dis- 
tribution. The two analysis methods may lead to very dif- 
ferent inferred correlation structure and interpretations. 

It is possible to apply the sibling-mixed model-fit pro- 
cedure to conditional angular correlations based on bin- 
ning (j> t ,pt) or (yt-,yt) as well. In that case the results 
may be very different from what is inferred by ZYAM 
subtraction, even when the pt cuts are identical. In par- 
ticular, jet yields may be much larger in the former case 
and the correlation shapes much different (27| . 



D. Relating the SS 2D peak to azimuth multipoles 

We have demonstrated that all azimuth multipoles can 
be predicted by measured properties of the SS 2D peak, 
AS dipole and nonjet quadrupole. The systematic vari- 
ation of just five model parameters controls all ID mul- 
tipole phenomenology in nuclear collisions, whatever the 
physical interpretation. "Higher harmonic flows" (trian- 
gular flow, etc.) are determined completely by the three- 
parameter SS 2D peak. The 2D peak and the claimed 
flow system are therefore equivalent within the ID pro- 
jection. However, the three-parameter SS 2D peak (com- 
bined with nonjet quadrupole and AS dipole) is a com- 
plete representation of all ID multipoles, whereas the ID 
multipole system cannot reconstruct the SS 2D peak and 
is therefore incomplete. Based on full 2D correlations we 
conclude that all structure assigned in ID to "higher har- 
monics" has a strong 2D dependence (large curvature on 
77a) not expected from flow mechanisms (Sec. IV Cl) . 



XII. SUMMARY 

Two notable features have been reported in 2D [on 
(r], (/))} angular correlations at RHIC: a same-side (SS) 
2D peak elongated on pseudorapidity r\ in more-central 
Au-Au collisions and an away-side (AS) ID peak on cf> 
seemingly distorted after subtraction of a combinatoric 
background. The apparent AS distortion appears as a 
double peak near 7r, interpreted by some as evidence for 
Mach cone formation in a dense medium. 

Recently, explanations for both apparent distortions 
have been proposed based on conjectured initial-state ge- 
ometry multipoles combined with collective expansion to 
produce a long-range (on r)) final-state quadrupole (in- 
creased elliptic flow) and sextupole (triangular flow) in 
the final state. Initial-state multipoles are said to repre- 
sent fluctuations in the A-A overlap geometry; final-state 
multipoles arc described as higher harmonic flows. 

In this article we introduce the measured systemat- 
ics of the SS 2D peak, AS dipole and nonjet azimuth 
quadrupole on A-A centrality. We show that the proper- 
ties of those three structures predict all azimuth multi- 
poles that can be derived from Fourier analysis of the ID 
projection onto azimuth of 2D angular correlations. We 
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show that some analysis methods relying on ID projec- 
tion can confuse Fourier components of a ID Gaussian 
with nominally independent other elements (e.g., flow 
and nonflow). 

We have reviewed the recent introduction of triangular 
flow and higher harmonic flows as a proposed explana- 
tion for elongation of the SS 2D peak. We show that 
the measured 2D peak properties do predict all triangu- 
lar flow data V3 and higher harmonics U4, etc., both at 
200 GeV, and with minor modification at LHC energies. 
We also demonstrate that Fourier analysis alone cannot 
describe 2D angular correlations. The ID Fourier model 
is falsified by 2D data. The 2D peak itself is then a more 
efficient and complete representation of 2D angular cor- 
relations than the separate azimuth multipoles derived 
from Fourier analysis of the ID projection. 

The surprising abundance of minijets in FS A-A corre- 
lation structure has presented a serious problem for pro- 
ponents of hydro-dominated RHIC collisions. Observed 
minijet systematics are consistent with negligible loss of 
minimum-bias scattered partons to thermalization even 
in central Au-Au collisions, leading to a recent trend to 
reinterpret nominal jet correlations as flow manifesta- 
tions. Based on Fourier analysis of ID projections onto 
azimuth collisions are said to be flow-dominated, with mi- 
nor additional contributions characterized as "nonflow." 
We have demonstrated that ID projection and Fourier 
coefficients alone provide an incomplete data description 
susceptible to incorrect inferences. 

Any comprehensive description of high-energy nuclear 
collisions must accommodate two-dimensional pQCD jet 
structure in p-p and peripheral A-A collisions and admit 
the possibility of modified jet formation in more-central 
A-A collisions. The full structure of the same-side 2D 
peak on 77 should be acknowledged as distinct from other 
correlation components. The centrality dependence of 
individual correlation components must be considered. 
The particle yield associated with the same-side peak 
(nominally jet fragments predicted by pQCD) should be 
described quantitatively. The same-side peak does have 
well-defined Fourier components when projected onto ID 
azimuth, but interpretation of any single Fourier term as 
representing flow can be questioned. 

This work was supported in part by the Office of Sci- 
ence of the U.S. DOE under grant DE-FG03-97ER41020. 



Appendix A: correlation models 

2D angular correlations are defined on difference- 
variable space (t]a,4>a) within angular acceptance 
(A77, A0). Accurate tja dependence is essential to distin- 
guish among data model elements and hadron produc- 
tion models. ID projections abandon key information 
and greatly reduce the ability of projected data to test 
theoretical models. Any structure projected onto peri- 
odic ID </>a can be represented by a Fourier series, but 
interpretation of the coefficients may be ambiguous. 



1. Model functions and 2D histograms 

2D correlation structure within some limited 77 ac- 
ceptance can be separated into ^-independent and 77- 
dependent components, as in Eq. (@|. The SS 2D peak 
centered at the origin always has a large curvature on 
r/A well-described by a Gaussian. The two sinusoids (AS 
dipole and nonjet quadrupole) have negligible curvature 
on 77A within the STAR TPC acceptance. The strong 
curvature on 77A of the SS 2D peak contradicts any de- 
scription relying solely on azimuth multipoles uniform 
on ?7a- Each term in Eq. (U) is associated with a cor- 
responding correlation model component on transverse 
momentum (pt,Pt), which may help to reveal the source 
mechanism and falsify incorrect models. 

In elementary high-energy hadron collisions we expect 
significant jet structure, including a narrow SS 2D peak 
at the origin and an AS ID peak (ridge) at 4>a = 7r rep- 
resenting momentum conservation for IS pairs of large- 
angle-scattered partons 0,[l7j- The 2D data model which 
arose from empirical study of angular correlations in p-p 
and peripheral A-A collisions includes three major ele- 
ments. The minimum-bias jet contribution accounts for 
two terms in Eq. ([3]) [HI, EH . In more-central A-A colli- 
sions an additional nonjet azimuth quadrupole term in- 
dependent of 77A emerges as a third term. 

2. Optimum data model 

Is there a unique 2D data model? If not, can a "best" 
of several models be identified? According to Ockham's 
razor the smallest number of model parameters is pre- 
ferred. All model elements should be necessary and the 
overall model should be sufficient. Necessity of a term in 
Eq. ^ can be tested by omitting it individually from the 
fit model. If the corresponding shape appears in the fit 
residuals the term is necessary to the data. If a combina- 
tion of necessary terms completely describes the data, 
with no significant structure in the fit residuals, that 
fit model is sufficient. A necessary and sufficient (N-S) 
model is statistically equivalent to the 2D data and may 
then be the "best" (possibly unique) model. Uniqueness 
may be challenged, but competing models must describe 
all aspects of the 2D data, including p t dependence, cen- 
trality dependence and collision-energy dependence, as 
well as or better than the N-S model that is challenged. 

3. ID projections and Fourier series 

2D angular correlations can be projected onto ID 77A 
or 0a- The projection onto 4>a is conventionally de- 
noted dihadron correlations with ZYAM subtraction [2lJ ■ 
ID projection is also implicit in nongraphical numerical 
methods (conventional V2 analysis) [35j . Any 2D data 
projected onto periodic ID azimuth can be represented 
by a single Fourier series, but such ID projections may 
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superpose multiple 2D model elements onto each Fourier 
coefficient (e.g., flow and nonflow) leading to confused 
interpretations. 

Replacing the SS 2D Gaussian by a ID Fourier se- 
ries increases the number of model parameters and de- 
grades the 2D data description, and the ID Fourier se- 
ries cannot describe the unprojected 2D angular correla- 
tions. Adding additional Fourier terms to the N-S model 
of Eq. (|4]) would degrade the capacity of the model to 
test hypotheses. More than one model element may then 
represent the same data feature leading to redundancy. 
Large excursions of fitted parameters may give the mis- 
leading impression that the original N-S model is an in- 
sufficient data description. 



Appendix B: Flow and Nonflow 

The terms "flow" and "nonflow" present an impor- 
tant definition problem for angular correlations at the 
RHIC. Flow in that context refers to "elliptic flow," a 
conjectured hydrodynamic phenomenon [351 ] . Conven- 
tional descriptions of elliptic flow assume that projected 
azimuth correlations are well described by a total az- 
imuth quadrupole t>2 oc cos(20a) whose dominant contri- 
bution is elliptic flow. Nonflow is assumed to be a sum of 
small contributions to the total quadrupole from extrane- 
ous structures (e.g., resonance correlations and possibly 
jets) which contribute small biases to V2 measurements. 

Definition problems arise from a breakdown of basic as- 
sumptions. Minimum-bias angular correlations are actu- 
ally dominated by a jet-like structure clearly apparent in 
2D correlations which has a large curvature on tja not ex- 
pected for hydro phenomena. The jet-like structure pro- 
jected onto azimuth has a substantial quadrupole compo- 
nent in its Fourier representation. A second quadrupole 
component with negligible curvature on ?7a distinct from 
jet-like structure is also apparent. The two quadrupole 
contributions may be better described as jet-related and 
nonjet quadrupoles. However, they also correspond sys- 
tematically to the terms nonflow and flow respectively. 

Several methods have been invoked in attempts to re- 
duce nonflow contributions to measured Vi. Higher cu- 
mulants such as i>2{4} were expected to reduce sensitiv- 
ity to few-particle correlations (such as jets) but may 
also be biased (negatively) by conjectured flow fluctua- 
tions. Combinations of V2 measurements have been used 
in attempts to estimate nonflow bias [35j . A favored 
method relies on 77 exclusion cuts. An exclusion zone 
|?7a| < 0(1) is expected to eliminate nonflow given the 
assumption that nonflow is always a "short-range" phe- 
nomenon. With the appearance of the ridge in triggered 
dihadron correlations and 77 elongation of the SS 2D peak 
in minimum-bias correlations the possibility of jet struc- 
ture (nonflow) outside the rj exclusion zone has emerged. 

In a recent paradigm shift maximum possible v m val- 
ues are sought using i> m {2} and higher multipoles are 
included: "This [extension to higher m] completely re- 



moves any assumption about the nature of non-flow cor- 
relations, except that [sic] is negligible at large relative 
pseudorapidity. As we have seen, this is a well-supported 
assumption" That strategy marks a reversal of 

previous attempts to discriminate "elliptic flow" (nonjet 
quadruple) from "nonflow" (jet- related quadrupole). It 
is assumed that any deviation from jet structure observed 
in p-p collisions, or an idealized (notional) jet definition, 
is a manifestation of flows. 



Appendix C: 2D multipole correlations 

We consider the 2D consequences of manipulating ID 
azimuth multipoles. Simulations are based on measured 
minimum-bias angular correlations from central Au-Au 
collisions 0, H[ . Curvature effects are emphasized by ex- 
trapolating the data parametrization to pair acceptance 
VA € [-4,4]. The STAR TPC acceptance extends to 
VA € [—3, 3] but is typically limited to [-2,2] for correla- 
tion analysis to reduce small artifacts at the acceptance 
edges. 



1. Multipole decomposition of the SS 2D peak 

Figure [18] (left panel) shows simulated data (nominal 
jet angular correlations) from central 200 GeV Au-Au 
collisions (y = 6, panel (d) in Fig. [5]). The nonjet 
quadrupole component Aq{2D} for that centrality is es- 
sentially zero. There are then two correlation compo- 
nents, the SS 2D peak and the AS dipole. Figure [TS] 
(right panel) shows the left panel with the AS dipole re- 
moved to isolate the SS 2D peak. 




FIG. 18: (Color online) Left: Simulated angular correla- 
tions from 200 GeV Au-Au collisions with centrality v — 6 
(b — 0) extrapolated to r\ £ [—2, 2] with nonjet components 
(quadrupole, ID Gaussian on ?7a) subtracted to reveal jet 
correlations. Right: The same histogram with the away-side 
dipole term subtracted to isolate the same-side 2D jet peak. 

Arguments in Refs. [l(| and others actually relate 
to the SS 2D peak, as we have shown in Sec. |X] The SS 
peak for minimum-bias angular correlations is consistent 
with a single 2D Gaussian in more-central Au-Au colli- 
sions. The SS peak retains a large curvature on r/^ for 
all collision conditions. Figure Q1J] shows four multipoles 
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(D,Q,S,0) from a 2D Fourier decomposition of the SS 
peak. Those structures correspond to the ID projections 
in Fig. [7] (right panel, dotted curves for m = 1...4). 




FIG. 19: (Color online) Same-side 2D multipole components 
A X {SS} of the SS 2D peak in Fig. \M (right panel): X 
dipole D, quadrupole Q, sextupole S and octupole O. 



2. 2D ZYAM subtraction 

We can use the 2D multipoles to illustrate the effect on 
2D angular correlations of ID ZYAM subtraction. Fig- 
ure [20] (left panel) shows 2D correlations from Fig. [18] 
(left panel) after subtraction of 75% of the SS 2D peak 
quadrupole V2{SS}, typical V2 oversubtraction from non- 
flow bias in ZYAM subtraction. The AS dipole is some- 
what distorted, but the SS 2D peak is still prominent. 

Figure I2TJ1 fright panel) shows the ID projection of the 
left panel as the bold histogram labeled ZYAM. The base- 
line estimate for those data according to the ZYAM prin- 
ciple (the "background" to be subtracted) is denoted by 
the dotted line. The actual jet structure is indicated by 
the dash-dotted and dashed curves for SS and AS peaks 
respectively relative to true zero. The structure of the 
ZYAM histogram is typical of what emerges from such 
analysis [271 ] . The apparent magnitude of the jet yield is 
greatly underestimated as a result of the ZYAM analysis 
but is still clearly evident in 2D angular correlations. 

In Fig. [5T] we simulate the suggestion in Ref. [l(| to ex- 
tend the ZYAM procedure by subtracting multipoles for 
in = 2, 3, 4. In the left panel we subtract ID multipoles, 
in the right panel the 2D multipoles from Fig. [19] The 
three SS multipoles combined are approximately equiv- 
alent to the SS peak minus the SS dipole, for either the 
ID or 2D case. As a consequence of ID subtraction (left 




FIG. 20: (Color online) Illustration of 2D ZYAM subtrac- 
tion with ID quadrupole (75% of ID «|{<SS'} quadrupole am- 
plitude) Left: 2D angular correlations resulting from ZYAM 
subtraction applied to the data in Fig. 1181 (left panel), Right: 
ID projection of the histogram in the left panel (bold solid 
histogram), the actual SS peak (dash-dotted curve) and AS 
peak (dashed curve) with their sum (light solid curve). 




FIG. 21: (Color online) 2D histograms with ID or 2D multi- 
poles for m = 2. ..4 subtracted. Left: 2D histogram in Fig. 1181 
after subtracting three SS ID multipoles. The AS dipole is 
canceled and a negative SS ID peak is superposed on the 
positive SS 2D peak. Right: The same 2D histogram after 
subtracting three SS 2D multipoles. The SS 2D peak is can- 
celed and a positive SS 2D dipole is superposed on the AS ID 
dipole. Both histograms project to flat ID histograms on <^>a • 

panel) the AS dipole is canceled by the SS dipole and 
the projected SS ID Gaussian, extrapolated across the 
r\ acceptance, is subtracted from the SS 2D peak. The 
highly-structured 2D histogram, when projected onto ID 
azimuth, is flat (uniform on azimuth) with value « 0.3. 

As a consequence of 2D subtraction (right panel) the 
SS 2D peak is canceled and the SS 2D dipole is added to 
the AS ID dipole. The highly-structured 2D histogram, 
when projected onto azimuth, is again fat with value 
« 0.3. For this demonstration the AS dipole amplitude 
Ad was adjusted to just cancel the projected SS 2D peak 
dipole (Ad F1GA2D), resulting in flat ID projections 
after subtraction. In real systems the AS dipole ampli- 
tude tends to be somewhat larger relative to the SS peak. 
The same analysis would then result in a small net AS 
dipole, interpreted by some as "global momentum con- 
servation." By subtracting the three multipoles all evi- 
dence of jet structure can be removed from projected ID 
azimuth correlations. But the jet structure is still clearly 
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evident in 2D correlations, albeit strongly distorted. 



Intermediate expressions are 



Appendix D: r\ exclusion cuts 

We seek an expression for factor G(a VA ;r]i,T]2) relat- 
ing the amplitude of a 2D Gaussian projected onto </>a 
with 77 exclusion cuts denoted by (771, 772) imposed to the 
amplitude of its ID Gaussian projection. The cut sys- 
tem consists of accepted 77 intervals [—772, — f]i] and [771, 772] 
symmetric about the origin. Defining A77 = 7/2 — 771 and 
s = \/2a riA the projection factor is given by 



with 
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With C + D = A77 2 the other terms are given by 



A = F(2 m , 771 +772) -2771^(2771,771 +772) (D4) 
B = 2772^(771 +772,2772) - F(r?i +772,2772), 



and G(cr I)A ; 771, 772) is determined. The limiting case 
??i -»?2 corresponds to G(a VA ; 771, 772) -> G(a VA ; A77) 
in Eq. (Hi]). 
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